Optimal path of diffusion over the saddle point and fusion of massive nuclei 
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Diffusion of a particle passing over the saddle point of a two-dimensional quadratic potential is 
studied via a set of coupled Langevin equations and the expression for the passing probability is 
obtained exactly. The passing probability is found to be strongly influenced by the off-diagonal 
components of inertia and friction tensors. If the system undergoes the optimal path to pass over 
the saddle point by taking an appropriate direction of initial velocity into account, which departs 
from the potential valley and has minimum dissipation, the passing probability should be enhanced. 
Application to fusion of massive nuclei, we show that there exists the optimal injecting choice for 
the deformable target and projectile nuclei, namely, the intermediate deformation between spherical 
and extremely deformed ones which enables the fusion probability to reach its maximum. 
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in 



i 

C 

O 

o 



> 

in 
o 

o 

oo 

o 



X 



I. INTRODUCTION 



The saddle point passage problem is of great interest 
in various fields of physics, such as collision of molec- 
ular systems, atomic clusters, biomolecules, and so on. 
The previous studies on this issue were mostly concen- 
trated on a simple diffusive dynamics with single degree 
of freedom, where the Langevin equation with constant 
coefficients can be easily solved in the case of a quadratic 
potential However, since many processes obviously 
involve more than one degree of freedom, which the one- 
dimensional (ID) model does not distinctly hold, high 
dimension at least two dimension (2D) is necessary. A 
case in point would be the fusion reaction of massive nu- 
clei, where the fusion is induced by diffusion 0, H[ and 
the asymmetrical or the neck degree of freedom of com- 
pound nuclei needs to be considered 0, Q. For these 
systems with the contact point of two colliding nuclei be- 
ing very close to the conditional saddle point, the poten- 
tial energy surface (PES) around the saddle point can be 
approximated to be a quadratic type. Under this approx- 
imation, Abe et al. || obtained an analytical expression 
for the multi-dimensional saddle-point passing probabil- 
ity. Some authors discussed quantum effect of the fu- 
sion probability by using the real time path integral [7| 
or the quantum transport equation [H, 1| , respectively. 
Boilley et al. studied the influence of initial distribu- 
tion upon the passing probability. Anomalous diffusion 
passing over the saddle point of ID quadratic potential 
was also discussed in Ref. [ll|. Nevertheless, it is not 
completely clear for the dynamical role of non-transport 
degrees of freedom. This might be very important for 
the quasi-fission mechanism in the fusion reaction, be- 
cause the average path of the fusing system in a multi- 
dimensional PES should be controlled by the off-diagonal 



components of inertia and friction tensors before the sys- 
tem arrives firstly at the conditional saddle point. 

Recently, theoretical calculations for the fusion bar- 
rier distribution, accounting for the surface curvature 
correction to the nuclear potential, have been presented 
by Hinde et al. [H [H, The geometrical effect 

significantly changes the near-barrier fusion cross sec- 
tion and the shape of the barrier distribution through 
an angle-dependent potential, where the target nucleus 
bears quadrupole and bexadecapole deformations while 
the projectile one is in a spherical shape. In these calcu- 
lations, the surface curvature correction to the sphere-to- 
sphere nuclear potential exerts influences upon the fusion 
probability through the height of fusion barrier. How- 
ever, the dynamical coupling effect of various deformative 
degrees of freedom needs to be added from the viewpoint 
of fusion by diffusion [3j. 

The primary purpose of this paper is to study the in- 
fluence of coupling between two degrees of freedom upon 
the passing probability. In Sec. II, we report the ana- 
lytical expression of the saddle-point passing probability 
by solving the 2D coupled Langevin equation with con- 
stant coefficients. In Sec. Ill, we discuss the effects of 
off-diagonal components of inertia, friction and potential- 
curvature tensors and then determine the optimal diffu- 
sive path. Sec. IV gives an application of this study to 
the actual fusion process of massive nuclei. A summary 
is written in Sec. V. 



II. THE PASSING PROBABILITY 



We consider the directional diffusion of a particle in a 



2D quadratic PES: U{x ll x 2 ) 



XiXj with i,j = 1,2 



and detuiij < 0, the motion of the particle is described 
by the Langevin equation 
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mijXj(t) + PijXj(t) + LO tj Xj(t) = £i(t) 



(1) 



2 



with Xj(0) = Xjo and Xj(0) — Vjo, where xio < 
and Vio > 0. Here and below the Einstein summation 
convention is used. The two components of the ran- 
dom force are assumed to be Gaussian white noises and 
their correlations obey the fluctuation-dissipation theo- 
rem (&(t)fj(*0) = k B Tmr*p k:j 5(t ~ t'), where k B is the 
Blotzmann constant and T the temperature. 

Assuming that xi-axis is the transport direction [uiu < 
0] , we write the reduced distribution function of the par- 
ticle for X\ while the variables x%(t), vi(t) and vz(t) are 
integrated out, 



W(x 1 ,t;x w ,X2o,Vio,v 20 ) 



1 



exp 



2tt<j Xi (t) 

(x 1 (t)-(x 1 (t))) 2 



2<(t) 



(2) 



Integrating over x\ from zero to infinity, we determine the 
passing probability over the saddle point \x\ = x 2 = 0] 

as 

/>oo 

P{t-x 10 ,X2o,vio,v 2 o) = / W(x 1 ,t;x 10 ,X2o,vio,v 2 o)dx 1 

Jo 



-erfc 



V2<J Xl (t) 



(3) 



Applying the Laplace transform technique to Eq. (1), 
we thus get x\(t) and its variance <J Xl (t) at any time, 



x x {t) = (i x (t)) 



2 rt 

E 

i=l 



Hiit - t%{t')dt', 



(4) 



t ft! 

dhHiit-tx) / dt2(Z i (t 1 )£ j {t2))H j (t-t 2 ), 

Jo 

(5) 

where the mean position of the particle along the trans- 
port direction is given by 



(xi(t)) = [Ci{t)x ia + G l+ 2(t)v M ] , 



(0) 



i=l 



which relates to the initial position and velocity. The 
time-dependent factors in Eq. (6) with exponential 
forms according to the residual theorem are Ci(t) — 
£ _1 [Ei(s)/P(s)] (i = 1 . . .4), the two response functions 
in Eqs. (4) and (5) read H x {t) = C- 1 [F 5 (s)/P(s)] and 
H 2 {t) = £- 1 [F 6 (s)/P(s)], where C' 1 denotes the inverse 
Laplace transform. The expressions of P(s) and Fi(s) 
(i = 1, ...,6) are written in the appendix. 



III. THE OPTIMAL DIFFUSIVE PATH 
A. The coupling effect of two degrees of freedom 

As one has known in the ID case, the passing proba- 
bility increases from to 1 when the initial velocity of 
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FIG. 1: The critical velocity (left) and the stationary passing 
probability (right) as functions of various off-diagonal compo- 
nents mi2, 0i2, and u>i2, respectively. The parameters used 
are: mu = 1.5, m22 = 2.0, /3u = 1.8, P22 = 1.2, W11 = —2.0, 
u>22 = 1.5, and 8 = 0. The initial velocities of the particle are 
v = 4.0, 2.2, 1.9, 1.6 from top to bottom (right). 



the particle increases. The critical velocity is defined by 
the passing probability being equal to \. This leads to 
the condition: limt_ >00 (a;i(t)) = 0. If all the off-diagonal 
components of the three coefficient tensors are not con- 
sidered, as well as X20 an d ^20 are taken to be zero, 
the critical velocity is determined from Eq. (6): Vq = 
[Fi(a)/F3(a)]xio, where a is the largest positive root of 
P(s) = 0. This is in fact identical to the one-dimensional 
result: w£ = -xio(v^ + iuj 11 +(3 11 )/(2m 11 ) [6], which 
is proportional to the friction strength. 

We now consider all the off-diagonal components of 
three coefficient tensors, namely, the correlations of two 
degrees of freedom are taken into account, the critical 
velocity can also be determined by ]ha.t—*oo{%i(t)) = 
and results in 



Cj (00)2:10 + C 2 (oo)x 2 o 
03(00) cos 9 + (74(00) sin( 



(7) 



where 9 denotes the incident angle between the initial 
velocity and the xi-direction, hence v w = v cos9 and 
V20 = v sin9. 
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FIG. 2: The stationary passing probability as a function of the 
off-diagonal component W12 for various /3i2. The parameters 
used are the same as those in the figure 1. 



FIG. 3: The stationary passing probability as a function of the 
incident angle. Here mn = 0.6, /3i2 = 0.8, and CJ12 = —0.5, 
as well as the other parameters are used as those in the figure 
1. 



Quantities plotted in the figures 1-7 are dimension- 
less and k B = 1.0 except for the units having been in- 
cluded in the figure caption. In Fig. 1, wc plot the 
critical velocity and the stationary passing probability 
as functions of the off-diagonal components of three co- 
efficient tensors, where one of off-diagonal components 
varies and the other two are fixed to be zero. The 
stationary passing probability is calculated by P paS s = 
linit-^oo icrfc [— (xi(i)) /(y/2a Xl (t))l . It is shown that the 
critical velocity increases with increasing the absolute 
value of mi2 or W12; while decreases with the increase 
of |/?i2|- The larger critical velocity of a system needs, 
the more difficult for the particle is to arrive at the po- 
tential top. This also implies that the passing probability 
is small when the dissipation along the diffusive path is 
large if the potential differences between the saddle point 
and the initial positions are equivalent. As is shown in 
the figure, the behavior of the passing probability is op- 
posite to that of the corresponding critical velocity. 

Figure 2 shows the stationary passing probability in 
the presence of two off-diagonal components U12 and /?i2, 
simultaneously, for mi2 = 0. It is seen that the maximum 
of the passing probability does not appear in the vertical 
case (u)\2 = 0). In the 2D PES, the particle is usually 
supposed to travel along the potential valley and then 
the steepest decedent direction. Because this is the di- 
rection which faces a smaller potential barrier. However, 
it may not be a path with a weaker damping. Under 
the effect of the off-diagonal component of the friction 
tensor, the particle is forced to select a better path with 
both low potential barrier and weak friction to surmount 
the saddle point of the potential. 



B. Determination of the optimal path 

Where is the optimal incident direction which enables 
the particle with given initial kinetic energy to have a 
larger passing probability? In order to determine this 
direction, we need to choose a special angle 9 m which 
enables the critical velocity to reach its minimum, i.e., 
from Eq. (7), 



dv o 

de 



0, 8 m — arctan 



C 4 (oo) 



(8) 



In fact, the largest analytically root of P(s) = domi- 
nantly determines the passing probability. The optimal 
incident angle 8 m can then be expressed by the Langevin 
coefficients as 

9 m = arctan {^^+^}-^^ + ^)\ . (9) 



m\\Fz{a) + TOi2F 6 (a) 

Using the same parameters as those written in Figs. 1 
and 3, we obtain 9 m ~ 0.258 rad, as is explicitly shown in 
Fig. 3, corresponding to the maximum of the stationary 
passing probability. 

We now define 7, ip and a to be the rotation angle 
of the major axis of the potential-curvature, friction and 
inertia tensors, respectively. They are found to have the 
following expressions: 



tan27 = 



tan2a 



2wi 



0^22 — Wn 

2mi2 
m 22 - mu ' 



tan2?/> 



2/3; 



12 



022 - Pll ' 



(10) 



As an example, for the case we have studied in Figs. 1 
and 3, these angles are: 7 ~ -7.973°, ip ~ -34.722°, 
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FIG. 4: (color online.) (a) The 2D potential energy sur- 
face, where the dotted curve is the saddle ridge line, (b) 
A schematic illustration of the optimal diffusion path (OP), 
where r m , r@, and r„ denote the major axes of the inertia, 
friction, and potential-curvature tensors, respectively. 



and a ~ 33.690°. For comparison, we write the optimal 
incident direction of the the particle we have obtained in 
the unit of one degree, i.e., 9 m ~ 14.779° (0.258 rad). 

In Fig. 4, we plot the two-dimensional quadratic po- 
tential and the optimal path in the x\-xi plane in a way of 
schematic illustration. All the coefficient elements used 
here have been written in the figures 1 and 3. It is illus- 
trated that the direction of the optimal path departs from 
the xi-direction. The effect of off-diagonal component of 
the inertia tensor makes the average path of the diffu- 
sive system turn toward the positive X2-axis, while the 
off-diagonal component of friction leads the mean path 
of the particle toward the negative a^-axis. Finally, the 
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FIG. 5: The stationary passing probability as a function of 
9 for various X20- Here Xio = — 1.0, vq — 1.9, mi2 = 0.6, 
P12 — 0.8, and u>i 2 = —0.5, as well as the other parameters 
are the same as the figure 1. 

competition of these two effects results in the optimal 
diffusive path shown in Fig. 4 (b). This phenomenon is 
similar to the quasi-stationary flow passing over the bar- 
rier in the fission case [l5[ , where the magnitude of the 
current is strongly influenced by the off-diagonal compo- 
nents of inertia and friction tensors. 

Figure 5 shows the dependence of the stationary pass- 
ing probability on the incident angle of the particle start- 
ing from various initial positions but with fixed initial 
kinetic energy. It is seen that the passing probability of 
the particle starting from a large positive X20 position is 
larger than that of starting from both small and nega- 
tive X20 positions. This is because the energy difference 
between the potential top and the initial position of the 
particle is small for the former. Amusingly, we find that 
the difference between the passing probabilities of two 
symmetrical positions (—1.0, —0.5) and (—1.0, 0.5) is ob- 
servably large. 

For a clearly understanding of the above results, we 
plot in Fig. 6 the mean diffusive path of a particle start- 
ing from different initial positions with different incident 
angles. Here (xi(t)) has been obtained in Eq. (6) and 

2 

(x 2 (t)) = [Ci+4(t)x i0 + C i+6 {t)v i0 ] , (11) 

i=l 

where all the time-dependent quantities are given in the 
appendix. The critical velocities are calculated by us- 
ing Eq. (7): = 1.5791 when x 20 = 0.5; v§ = 2.2321 
when X20 = —0.5, for x±q = —1.0 and 8 = 0. Hence 
the stationary passing probability of the particle starting 
from X20 = 0.5 is larger than that of the particle start- 
ing from X20 — — 0.5. In particular, under the present 
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FIG. 6: (color online.) The mean diffusion path of a particle 
starting from the initial positions of two kinds at fixed Vq = 
1.9, where the solid and dashed lines correspond to 8 = and 
6 = 0.258 rad, respectively. Here all the Langevin parameters 
are the same as the figures 1 and 5. 



circumstance, the diffusive process of the particle with 
different incident angles shows an interesting behavior. 
Because the initial velocity of the particle along the re- 
direction for 8 = is larger than that of 8 — 0.258 rad, 
the former can move to a position being closer to the 
saddle point than the latter. Thus the passing probabil- 
ity for 8 = is larger than that of 8 = 0.258 rad at the 
beginning. However, the width of the Gaussian distribu- 
tion is independent of the incident angle and increases 
with the increase of time. Although the center position 
of the particle's distribution with 8 = 0.258 rad is behind 
that of 8 = 0, as the time goes on, it will have a lager 
share of distribution passed the saddle point. Therefore, 
the passing probability for a particle with incident angle 
8 = 0.258 rad is larger than that of 8 = in the long 
time. 

Time-dependent passing probabilities shown in Figs. 
7 (a) and 7 (b) are also in complete agreement with the 
above theoretical analysis. 
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FIG. 7: (color online.) Time-dependent passing probability 
for various X20 an d 9. Here xio = —1.0, vo = 1.9, X20 = —0.5 
in (a) and £20 = 0.5 in (b), as well as the Langevin parameters 
are the same as the figures 1 and 5. 



IV. APPLICATION TO FUSION OF MASSIVE 
NUCLEI 

We now apply the present 2D simplified diffusive model 
to investigate the fusion of two massive nuclei, which has 
been described by directional diffusion over the saddle 
point @. As a particular example, we calculate the fu- 
sion probability of the nearly symmetrical reaction sys- 
tem 100 Mo+ 110 Pd 16], which is plotted as a function of 
the center-of-mass energy E c . m . in Fig. 8. A schematic il- 
lustration of the deformation of the compound nucleus is 
also shown in this figure. The temperature of the fusing 
system is determined by aT 2 — E CAn . + Q — Eb, where 
a = A/10 is the energy level constant with A the nucleon 
number of the compound nucleus, Q denotes the reaction 
Q value and Eb the barrier height of fission potential. 

The {c, h, a} shape parametrization [17] with the elon- 
gation c (the half of the nuclear length) and neck variable 
h are used, i,e, x\ = c, x% = h, and the asymmetrical pa- 
rameter a is fixed to be zero. The inertia and friction ten- 
sors are calculated by the Werner- Wheeler method and 
the one-body dissipative mechanism [l8[ , respectively, all 
the Langevin coefficients are considered to be constants 
at the saddle point. The three components of potential- 
curvature tensor are: u>n = —28.2304, W22 = 275.4211, 
and ui±2 — 50.4551 in the unit of MeV; the compo- 
nents of friction tensor /3 U = 701.9967, /3 22 = 621.4425, 
012 = 601.4934 in the unit of 10~ 21 MeV-sec; the in- 
ertia elements ran — 102.4081, 77122 = 134.4673, and 
TO12 = 110.3783 in the unit of 10~ 42 MeV-sec 2 . 

It is found a highlighted interesting result from Fig. 8 
that there exists the optimal collision shape for projec- 
tile and target nuclei, which induces the maximum fusion 
probability under the same center-of-mass energy. It can 
be easily understood from a combining role of the off- 
diagonal components of three dynamical coefficient ten- 
sors. This concludes that the fusion probability of mas- 
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FIG. 8: The fusion probability of the reaction 100 Mo + 110 Pd 
as a function of the center-of-mass energy for various initial 
positions. Schematically illustrated as well are the deformed 
shape of the compound nucleus. 



FIG. 9: The fusion probability of the reaction 100 Mo + 110 Pd 
as a function of the center-of-mass energy for the situations 
with and without off-diagonal components. 



sive nuclei can be enhanced if the two collision heavy ions 
are polarized to be ellipsoid and the collision direction be- 
tween the long and short axis of ellipsoid is appropriately 
selected. For the fusion of deformed massive nuclei, there 
exists the optimal angle for the incident nucleus to col- 
lide with the target one, which favors the fusion of heavy 
ions to be accomplished. 

Figure 9 shows the fusion probability of 100 Mo+ 110 Pd 
as a function of the center-of-mass energy when the off- 
diagonal components are considered partly. Here the ini- 
tial position of the fusing system is chosen into the opti- 
mal one, i.e., cq — 1.75 and ho = 0.2789 as well as all the 
parameters used are the same as the figure 8. This fusion 
reaction system is regarded as an example to compare the 
results with and without off-diagonal components in the 
potential surface, friction and mass parameters, which 
is reflected in the result presented in the above section. 
In fact, the case of without all off-diagonal components 
is equivalent to the one dimension or without the neck 
situation 

It is seen from Fig. 9 that the increase of the ID fusion 
probability curve versus the energy is faster than that 
of 2D case. It has been known that the pervious ID 
diffusion model without the neck variable proposed the 
fusion probability larger than the experimental data, so 
the present completely coupled 2D diffusion model might 
be appropriate. Moreover, the results for the presence 
of only one of three off-diagonal components can also 
be understood from the critical velocity (kinetic energy) , 
to see Fig. 1. Namely, the larger the critical kinetic 
energy of the system is, the less the fusion probability 
is for the same center-of-mass energy. The nonvanishing 
/3i2 allows the smallest critical kinetic energy and the 



presence of mi 2 leads to the largest critical kinetic energy. 
Therefore, we have a relation for the fusion probabilities: 
Pf us (/3i2 ^ 0) > P fus (wi2 ^ 0) > Pf us (mi2 == 0) at a fixed 
center-of-mass energy. 



V. SUMMARY 

We have studied the diffusion process of a particle 
passing over the saddle point of a two-dimensional non- 
orthogonal quadratic potential. The expression of the 
passing probability is obtained analytically, where the in- 
ertia and friction tensors are not diagonal. The optimal 
incident angle of the particle's initial velocity is deter- 
mined. Our results have shown that the optimal diffusive 
path, which departs from the potential valley in the two- 
dimensional potential energy surface, induces the maxi- 
mum saddle-point passing probability. This is due to the 
competition effect between the off-diagonal components 
of inertia, friction and potential-curvature tensors. Wc 
have investigated the fusion probability of massive nuclei 
and compared the results with and without off-diagonal 
terms, for instance, the reaction of 100 Mo-|- 110 Pd. Due to 
the influences of off-diagonal components of inertia and 
friction upon the diffusive path, which are calculated by 
the {c, h, a} parametrization, the fusion probability can 
be enhanced for an appropriate choice of the collision 
direction of the deformable target and projectile nuclei. 
The optimal configurations of colliding nuclei is between 
spherical and extremely deformed ones. In further, the 
present study also provides useful information in connec- 
tion with the synthesis of superheavy elements. 
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The quantities appear in the expression of {x\{t)) are 

P(s) = (detm)s + (mnfe + m 22 /3n - 2m 12 (3i 2 )s A 

+(det/3 + miiLu 22 + m 22 wn — 2mi 2 wi 2 )s 2 

+(/3uw 22 + /3 22 wn - 2/3i 2 wi 2 )s + detw, 
Fi(s) = (detm)s 3 + (mnfe + m 22 /3n - 2m 12 f3 12 )s 2 

+(det/3 + m n uj 22 - m 12 u 12 )s + /3nw 22 - [3 12 uj 12 , p 10 (s) 
F 2 (s) = (mi 2 uj 22 - m 22 uji 2 )s + (3i 2 uj 22 - f3 22 LUi 2 , 
F 3 (s) = (detm)s 2 + (mnfe — m\ 2 {5\ 2 )s + mnLu 22 

-mi 2 wi 2 , 

F 4 (s) = {m 12 (3 22 ~ m 22 (3 12 )s + m 12 u 22 - m 22 uj 12 , 
F 5 (s) = m 22 s 2 + (3 22 s + uj 22 , 

F 6 (s) = -m 12 s 2 - f3 12 s-u> 12 , (12) 



